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We calculate the probability ("quenching weight") that a hard parton radiates an 
additional energy fraction AE due to scattering in spatially extended QCD mat- 
ter. This study is based on an exact treatment of finite in-medium path length, it 
includes the case of a dynamically expanding medium, and it extends to the an- 
gular dependence of the medium-induced gluon radiation pattern. All calculations 
are done in the multiple soft scattering approximation (Baier-Dokshitzer-Mueller- 
Peigne-Schiff-Zakharov "BDMPS-Z" -formalism) and in the single hard scattering 
approximation (N = 1 opacity approximation). By comparison, we establish a sim- 
ple relation between transport coefficient, Debye screening mass and opacity, for 
which both approximations lead to comparable results. Together with this paper, a 
CPU-inexpensive numerical subroutine for calculating quenching weights is provided 
electronically. To illustrate its applications, we discuss the suppression of hadronic 
transverse momentum spectra in nucleus- nucleus collisions. Remarkably, the kine- 
matic constraint resulting from finite in-medium path length reduces significantly 
the pj_-dependence of the nuclear modification factor, thus leading to consistency 
with the data measured at the Relativistic Heavy Ion Collider (RHIC). 



1. Introduction 



Hard partons produced in nucleus-nucleus collisions at RHIC and LHC prop- 
agate through highly excited matter before hadronizing in the vacuum. The 
resulting medium-dependence of parton fragmentation is expected to affect 
hadronic observables. This is of two- fold interest. First, it provides a novel 
test of the space-time evolution of the perturbative parton shower. Second, 
the modification of hadronic observables due to the spatially extended, hot 
and dense QCD matter allows to characterize the properties of the transient 
state produced in the collision. 

Gluon emission off highly virtual hard partons is an essential component in 
the standard description of parton fragmentation in elementary processes. This 
effect degrades the energy of the leading parton. Recently, it has been pro- 
posed 1 that in the presence of a spatially extended medium, the additional 
medium-induced energy degradation of the leading parton can be described 
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by a probability P(AE), the so-called quenching weight, which is obtained 
from a probabilistic iteration of the medium-modified elementary splitting 
processes q — > qg and g — > gg. The main purpose of the present work is to cal- 
culate and compare this quenching weight for different approximations of the 
medium- modified splitting process, to make the results for P(AE) available 
as a numerical subroutine, and to illustrate the use of this subroutine with 
some applications. 

We start from recent calculations 2 5 of the modification of the elementary 
splitting processes q — > qg and g — > gg due to multiple scattering. These 
results go under the name medium-induced gluon radiation. They present 
limiting cases of an unique path-integral expression given in eq. (2.1) below. 
Technically, they collect all terms to leading order in nuclear enhanced modi- 
fications 0(a s A 1 ^ 3 ), thus accounting for the leading additional interactions of 
the parton shower with the medium. 

The paper is organized as follows: In section 2, we compare the medium- 
induced gluon energy distribution radiated off a hard parton in two limits 
which emphasize the role of multiple soft and single hard medium-induced 
scatterings, respectively. In section 3 we give results for the quenching weights 
corresponding to these limits. These quenching weights can be calculated with 
the numerical subroutine accompanying this paper. In section 4, we extend 
these calculations to the case of an expanding medium, and in section 5, we 
discuss the extension to radiation within a finite cone. As application, we 
calculate in section 6 in two different approaches the suppression of hadronic 
transverse momentum spectra and we compare our results to the nuclear mod- 
ification factor measured in Au-Au collisions at the Relativistic Heavy Ion 
Collider (RHIC). 



2. Medium-induced gluon radiation from a static medium 



The inclusive energy distribution of gluon radiation off an in-medium produced 
parton takes the form 4 ' 6 ' 7 
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(2.1) 



Here, denotes the transverse momentum of the emitted gluon. The limit 
k± = |kj_| < x°° 011 the transverse phase space allows to discuss gluon emis- 
sion into a finite opening angle 6, x — sin 6. For the full angular integrated 
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quantity, x — 1- 



The two-dimensional transverse coordinates u, y and r emerge in the deriva- 
tion of (2.1) as distances between the positions of projectile components in the 
amplitude and complex conjugate amplitude. The longitudinal coordinates yi, 
jji integrate over the ordered longitudinal gluon emission points in amplitude 
and complex conjugate amplitude, which emerge in time-ordered perturba- 
tion theory. These internal integration variables play no role in the following 
discussion. They are explained in more detail in Ref. 4 . 

The radiation of hard quarks or gluons differs by the Casimir factor Cr — Cf 
or Ca, respectively. Numerical results are for fixed coupling constant a s = 1/3, 
except where stated otherwise. The properties of the medium enter eq. (2.1) in 
terms of the product of the time-dependent density n(£) of scattering centers 
times the strength of a single elastic scattering er(r). This dipole cross section 
er(r) is given in terms of the elastic high-energy cross section |a(q)| 2 of a single 
scatterer, 



In this section, we study the energy distribution (2.1) for a static medium in 
the limiting cases of multiple soft and single hard momentum transfer. Then 
we compare these two limiting cases of (2.1). 

A. Multiple soft scattering approximation 

For arbitrary many soft scattering centers, the projectile performs a Brownian 
motion in transverse momentum. This dynamical limiting case can be studied 
in the saddle point approximation of the path-integral (2.1), using 3 ' 8 



Here, q(£) is the transport coefficient 9 which characterizes the medium- 
induced transverse momentum squared {qj_) mc d transferred to the projectile 
per unit path length A. For a static medium, the transport coefficient is time- 
independent, 
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In the approximation (2.3), the path integral in (2.1) is equivalent to that of a 
harmonic oscillator. The corresponding analytical expressions are summarized 
in Appendix A. 

Qualitative arguments 10 : We consider a gluon in the hard parton wave func- 
tion. This gluon is emitted due to multiple scattering if it picks up sufficient 
transverse momentum to decohere from the partonic projectile. For this, the 
average phase f accumulated by the gluon should be of order one, 

V =/£a*\~^L=^. (2.5) 
r \2uj I 2uj uj K J 



Thus, for a hard parton traversing a finite path length L in the medium, the 
scale of the radiated energy distribution is set by the "characteristic gluon 
frequency" 

u c = ^qL 2 . (2.6) 



For an estimate of the shape of the energy distribution, we consider the number 
iV co h of scattering centers which add coherently in the gluon phase (2.5), k\ ~ 
^coh (<?l)mcd- Based on expressions for the coherence time of the emitted gluon, 
t coh ~ pj- ~ and N coh = ^ = y ^^V^ , one estimates for the gluon 

energy spectrum per unit path length 

dl 1 dl lscatt a s [J 

u ~y — t — iTr — u ~y — ; — — — a s \ — ■ (2-7) 

dcu dz N coh du) dz t co h V w 



This 1 / ^/cJ-energy dependence of the medium-induced non-abelian gluon en- 
ergy spectrum is expected for sufficiently small uj < uj c . 

Quantitative analysis: The gluon energy distribution (2.1) depends not only 
on uj c , but also on the constraint kj_ < on the transverse momentum phase 
space of the emitted gluon. This enters the calculation via the dimensionless 
kinematic constraint 11 

R x = X -q X 2 L\ R = R x=1 = iu c L. (2.8) 



This constraint is neglected in the argument leading to the 1 / ^/uJ-energy de- 
pendence of (2.7). In the following sections, we limit the discussion to angular 
fully integrated quantities for which \ — \. The only exception will be the 
discussion of the angular O-dependence of in section 5, where we use 
X = sin 0. 
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The limit R — > oo which removes the kinematic constraint from (2.1) is ei- 
ther realized by extending the /^-integration ad hoc to infinity. Alternatively, 
R — > oo can be viewed as the limit of infinite in-medium path length since 
it corresponds to L — > oo for x an d uo c fixed. In Appendix A, we derive the 
R — > oo limit of the energy distribution (2.1), 



hm uj— = m 



COS 



(2.9) 



This coincides with the result of Baier, Dokshitzer, Mueller, Peigne and Schiff 9 . 
As expected from the estimates in (2.5) and (2.7), it shows a characteristic 
1/ Y^-energy dependence for small uj which is suppressed above the charac- 
teristic gluon frequency uj c -} 



.. dl 2a s C R J for UJ<UJ c, 

hm cj— ~ < / X 2 (2.10) 

[ _L (^j for > cu c . 



The average parton energy loss is the zeroth moment of this energy distribu- 
tion, 

(AE)^ = lim duu^- = ^u c . (2.11) 



This is the well-known /^-dependence of the average energy loss 2 ' 9,12 . Due to 
the steep fall-off at large u, the a;-integral in (2.11) is dominated by the region 
uj < uj c j \/2. 

We have evaluated numerically the energy distribution (2.1) for finite values 
of the density parameter R. As seen in Fig. 1, the distribution approaches for 
any value of R the BDMPS limit (2.9) at sufficiently large gluon energy. Below 
a critical gluon energy uj, however, the finite size gluon spectrum is depleted 
in comparison to the BDMPS limit. To understand this effect, we consider 
the characteristic angle 6 C at which medium-induced gluons are radiated on 
average, 

2 _ . k\ _ y/uq _ ( uj \ ~ 3 / 2 1 



UJ 2 UJ 2 \UJrJ R 



For C ~ 1, the emitted gluons are sensitive to the kinematic constraint since 
k± ~ 0(uj). The condition @ c ~ 1 thus provides an estimate for the gluon 
energy uj below which the energy distribution is cut off, 

13 « (If . ,2,3) 



uj c \RJ 
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The position of the maximum of uoj^ as a function of R is consistent with 
this dependence on lo, see Fig. 1. In general, gluon radiation at small energies 
corresponds to gluon radiation at large angle and is depleted as soon as the 
finite size of the transverse momentum phase space becomes relevant. This 
suppression of the non-perturbative small-u; contributions helps to make the 
calculation of medium-induced energy loss perturbatively stable. 



. BDMPS (R = oo) 




Fig. 1. The medium-induced gluon energy distribution 

^Tijjj * n ^6 multiple soft 

scattering approximation for different values of the kinematic constraint R = uj c L. 



The gluon energy distribution (2.1) also determines the multiplicity N(uj) of 
gluons emitted with energies larger than uj 



N(u) = f 

J Ul 



, dl(u') 
du> 



(2.14) 



In the absence of kinematic constraints, and for sufficiently small energies uj, 
one finds from (2.10) that the total multiplicity diverges like 

hm iV(ccO = ^^^, for u;<u c . (2.15) 

However, realistic kinematic constraints on the transverse momentum phase 
space (R < oo) deplete the gluon energy distribution at small uj and ensure 
that the total gluon multiplicity N(u = 0) is finite, see Fig. 2. 
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For realistic kinematic constraints, R < 10000, the average additional total 
multiplicity is N(u = 0) < 3. In comparison to the typically ~ 5 — 10 semi 
hard partons which are the partonic final state of a 100 GeV jet simulated in 
a parton shower, this additional multiplicity is not negligible. It supports the 
naive expectation that the number of partons in the jet increases and softens 
with increasing transport coefficient or path length. 




Fig. 2. The multiplicity of additional medium-induced gluons (2.14) radiated 
with energy large than uj. Calculation done in the multiple soft scattering approxi- 
mation. 

B. Single hard scattering approximation 

In the previous subsection, we have studied the energy distribution (2.1) of 
medium-induced gluon radiation in the limit in which the partonic projectile 
performs a transverse Brownian motion due to multiple soft scattering. Now, 
we consider the opposite limiting case in which the radiation pattern results 
from an incoherent superposition of very few tiqL single hard scattering pro- 
cesses positioned within path length L. This limit is obtained by expanding 
the integrand of the energy distribution (2.1) in powers of (n(£) a(r)) N up to 
first order 4 ' 5 ' 13 . Analytical expressions are given in Appendix B. 

Qualitative arguments: We consider a hard partonic projectile which picks up 
a single transverse momentum fi by interacting with a single hard scatterer. 
An additional gluon of energy uj decoheres from the projectile wave function 
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if its typical formation time t co h = j^r is smaller than the typical distance L 
between the production point of the parton and the position of the scatterer. 
The relevant phase is 

^ ^coh ^ ^ 



which indicates a suppression of gluons with energy to larger than the charac- 
teristic gluon energy 

Gj c = X -^L. (2.17) 



The gluon energy spectrum per unit path length can be estimated in terms of 
the coherence time t co h and of the average number n L of scattering centers 
contributing incoherently 

dI N=l a u 2 

u—— ~ (n L) p- ~ (n L)a s — . (2. 18) 

duj dz t coh to 



This is the typical 1/cj-dependence of the non-abelian gluon radiation spec- 
trum in the absence of LPM-type destructive interference effects. 

Quantitative analysis: We have calculated the first order in opacity n L of the 
gluon energy distribution (2.1). To first order, the entire medium-dependence 
comes from the interaction of the hard parton with a single static scatter- 
ing center, multiplied by the number n L — L/X of scattering centers along 
the path. Modeling the single scatterer by a Yukawa potential with Debye 
screening mass /i, we derive in Appendix B 

dI N=1 _ n a s C R t rA f°° r - sin(r) 



duo n Jo r 

x — , 1 : 1 ■ (2.19) 

V r + ^ J((R/2-y) +r + 1 ) 2 -4rR/2 1/ 



This energy distribution depends on the phase factor 7 defined in (2.16), and 
on the kinematic constraint in transverse momentum phase space, 

R x = \ xV L 2 , R = R x=1 = uj c L . (2.20) 

In what follows, we work for x — 1 except for the discussion of the angular 
dependence in section 5. In the limit in which the kinematic constraint is 
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removed, the characteristic l/cu-energy dependence of the estimate (2.18) is 
recovered for sufficiently large gluon energies uj > uj c , 
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(2.21) 



This limit agrees with the results of Gyulassy, Levai and Vitev 13 . The average 
parton energy loss for a single hard scattering is dominated by contributions 
from the region uj > uj c , 13,14 

dujuj—, ~ ^^(n L)uj c \og\E/uj c ] . (2.22) 

duj 2 

It is logarithmically enhanced in comparison to the region uj < uj c for which 
r r c A dlN=1 2a sC R f 

Jim / dujuj — ~ (noL)uj c . (2.23) 

R^oo JO duj 7T 
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Fig. 3. The medium-induced gluon energy distribution ujj^ for a hard quark 
in the single hard scattering approximation, calculated for different values of the 
kinematic constraint R. 
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Remarkably, the average parton energy loss receives its dominant contribution 
from the hard region u > lj c in the opacity approximation (2.22) but from the 
soft region u < oj c in the multiple soft scattering approximation (2.11). 

We have evaluated numerically the energy distribution (2.19) for finite values 
of the kinematic constraint R. In close analogy to the multiple soft scatter- 
ing approximation, the emission of soft gluons is suppressed in the opacity 
approximation due to the kinematic constraint R = uj c L on the transverse 
momentum phase space, see Fig. 3. To estimate the scale Co at which this 
suppression sets in, we parallel the argument leading to eq. (2.12). We require 
that the characteristic angle of the gluon emission is of order one, finding 




(2.24) 



The numerical position of the maximum of u duj in Fig. 3 changes oc in 
accordance with this estimate. We thus have a semi-quantitative understand- 
ing of how phase space constraints deplete the non-perturbative soft region of 
the medium-induced gluon energy distribution. 




C0/(a) c 

Fig. 4. The multiplicity of additional medium-induced gluons (2.14) radiated 
with energy larger than uj. Calculation done in the single hard scattering approxi- 
mation. 



In Fig. 4, the additional medium-induced gluon multiplicity (2.14) is calculated 
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in the opacity approximation. In the absence of kinematic constraints R = oo 
(R = oo) and for sufficiently small gluon energies u> < lj c (uj < cu c ), this 
multiplicity changes oc ^ in the opacity approximation (oc in the multiple 
soft scattering approximation). In the presence of kinematic constraints, the 
total additional multiplicity is comparable for both approximations: N(u = 
0) < 3. 



C. Comparison: multiple soft vs. single hard scattering approximation 

Qualitative: The transverse momentum accumulated by a projectile due to 
Brownian motion increases linearly oc ji 2 n^L with path length where no = j 
denotes the longitudinal density of scattering centers. This leads to 

fj, 2 n L = qL , for Brownian motion, (2.25) 
and thus 

to c = -qL 2 = — cu c , in the multiple soft scattering limit. (2.26) 
Z A 

Recent applications of the opacity approximation use 1 < j < 3. In this case, 
the gluon energy distribution is much harder in the opacity approximation 
than in the multiple soft scattering approximation, see Fig. 5. 

Quantitative: The relation u c = (n L) cu c holds only if the projectile accumu- 
lates transverse momentum by Brownian motion. In general, deviations from 
Brownian motion are due to the high transverse momentum tails of the elastic 
scattering cross sections, 

W^SGPW (2 ' 27) 

In QED, the Coulomb scattering distribution is well represented by the theory 
of Moliere and shows logarithmic deviations from Brownian motion. For QCD, 
one can identify an analogous logarithmic term in the transport coefficient 
(2.3) by expanding the dipole cross section (2.2), 



^ L=n " L \w? Kq)|2 \ q2 cos2 v ~ (n ° L) ^ 2 ln v if ■ (2 - 28) 



Here, E cut denotes the upper cut-off of the logarithmically divergent ^-integral. 
This changes eq. (2.26) to 
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(n L) u c In 
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(2.29) 



The logarithmic term makes the comparison between single hard and multiple 
soft scattering approximation more difficult. Based on (2.29), the curves for 
the single hard scattering approximation should be shifted in Fig. 5 by a 
factor In > 1 to the left. For realistic values \pi > Aqcd and E cut < E 



say], we find In \J^f- «C 10. Thus, the above conclusion stays unchanged: the 
medium-induced gluon energy distribution is significantly harder in the single 
hard scattering approximation than in the multiple soft one. 



BDMPS 

L/A=0.5 

L/A= 1 

L/A=2 

, L A = 4 




U)/ 0J C 



Fig. 5. The gluon energy distribution without kinematic constraint (R, R — > oo) 
as calculated in the multiple soft scattering approximation, and in the single hard 
scattering approximation for n$L = 0.5, 1, 2,4. Results for the single hard scattering 
approximation are plotted for (L/X)lo c = oj c . 



3. Quenching Weights 



Medium-induced gluon radiation modifies the correspondence between the ini- 
tial parton and the final hadron momenta. This modification can be deter- 
mined from the distribution P(AE) of the additional medium-induced energy 
loss which we calculate now. If gluons are emitted independently, P(AE) is 
the normalized sum of the emission probabilities for an arbitrary number of n 
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gluons which carry away the total energy AE: 
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• (3-1) 



The summation over arbitrarily many gluon emissions in (3.1) can be per- 
formed by Laplace transformation 1 



V{y) = exp 



L 



"*° , dliuj) 

dliJ : 

du 



(l - e-" w ) 



(3.2) 
(3.3) 



Here, the contour C runs along the imaginary axis with Rei/ = 0. In general, 
the probability distribution P(AE) has a discrete and a continuous part, 11 



P(A£) = poS(AE) +p(AE) . 



(3.4) 



The discrete weight po may be viewed as the probability that no additional 
gluon is emitted due to in-medium scattering and hence no medium-induced 
energy loss occurs. This weight is determined by the total gluon multiplicity, 



po = lim V{v) = exp [-N(u = 0)] . 



(3.5) 



For finite in-medium path length, there is always a finite probability po ^ 
that the projectile is not affected by the medium. Only a finite number of 
additional medium-induced gluons can be emitted, see (3.5). For infinite in- 
medium path length, one finds 

lim p = . (3.6) 

R— >oo 



The medium-induced gluon energy distribution uj^ determines to what extent 
the total energy distribution of a given parton deviates from its "vacuum" 
fragmentation in an elementary collision, 

d j(tot) rf/ (vac) d j 
OJ : — = UJ h UJ— . (3.7) 

duo duj duj 



From the Laplace transform (3.2), we obtain for the corresponding total prob- 
ability 

p( tot \AE)= dE P(AE - E) P^ C \E) . (3.8) 
Jo 
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The probability P( tot \AE) is normalized to unity and it is positive definite. 
In contrast, the medium-induced modification of this probability P(AE), is 
a generalized probability. It can take negative values for some range in AE, 
as long as its normalization is unity 

/ dEP(E)=p + dEp(E) = l. (3.9) 
Jo Jo 



In this section, we calculate P(AE) in the multiple soft and single hard scat- 
tering approximations. The results of these quenching weights are available as 
a FORTRAN routine 15 . 



A. Quenching weights in the multiple soft scattering approximation 



By numerical evaluation of the Laplace transform (3.2), (3.3), we have calcu- 
lated quenching weights P(AE) for the medium-induced energy distribution 
(2.1) in the multiple soft scattering approximation. To motivate the range of 
parameter values studied in what follows, we relate the transport coefficient q 
to the in-medium path length L and the saturation scale Q s , w 

Q 2 S - QL. (3.10) 
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Estimates for the saturation scale are very uncertain but Q 2 S < (3 GeV) 2 may 
be considered as an upper bound at LHC. To discuss in-medium path lengths 
L up to twice a nuclear Pb radius, we thus have to explore the parameter 
space up to R < 40000. We choose a very small lower value, R = x 2uJ cL 3 = 1, 
in order to tabulate quenching weights for the radiation outside very small 
opening angles x — sin 6. All results will be given for energies in units of uj c . 
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Fig. 7. The continuous part of the quenching weight (3.4), calculated in the 
multiple soft scattering approximation for a hard quark (upper row) or hard gluon 
(lower row). 



The discrete weight po of the probability distribution P(AE) is plotted in 
Fig. 6 as a function of the kinematic cut-off R = u c L. It approaches unity in 
the absence of a medium (R — > 0) and it vanishes in the limit of infinite in- 
medium path length, see (3.6). Remarkably, p exceeds unity for small values 
R < 100. This indicates that there is a phase space region at very small 
transverse momentum, into which less gluons are emitted in the medium than 
in the vacuum. The "vacuum" gluon radiation is shifted to larger transverse 
momentum in the presence of a medium 6 . The decrease of the discrete weight 
for large R and its growth above unity for sufficiently small R both depend on 
the strength of the interaction between partonic projectile and medium. They 
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are thus more pronounced for gluons than for quarks. 



The continuous part p(AE) of the probability distribution (3.4) is shown in 
Fig. 7 as a function of the dimensionless energy fraction AE/uj c for different 
values of the kinematic constraint R. Increasing the density of the medium 
(i.e. increasing the transport coefficient q) or increasing the in-medium path 
length L corresponds to an increase of u c and R. The Figs. 6 and 7 specify 
how the probability that the parton loses an energy fraction AE changes 
with these medium properties. As expected from the normalization (3.9), the 
continuous part p(AE) shows predominantly negative contributions for small 
values R < 100 where the discrete weight p exceeds unity. 



In the limit R — > oo, the quenching weight was found to be fit very well by a 
two-parameter log-normal distribution 16 . This is a heuristic observation which 
is difficult to connect to the analytical structure of the gluon energy distribu- 
tion. Analytically, an estimate of the quenching weight can be obtained 1 in 
the limit R — > oo from the small-u; approximation uj^ oc ^= in eq. (2.10), 



papprox/ \ 
r BDMS \ e ) 



a 



7r a 
e 



where a = Tr-^-uJr 



7T Z 



(3.11) 



This approximation is known to capture 11 the rough shape of the probability 
distribution for large system size, but it has an unphysical large e-tail with 
infinite first moment / dee P^^si 6 )- Also, its maximum e max = ^p- grows 
stronger with the effective coupling a s Cr than the numerical result in Fig. 7. 
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Fig. 8. The discrete part po of the quenching weight calculated in the single 
hard scattering approximation for opacity uqL = 1. 
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B. Quenching Weights in the opacity approximation 



We have evaluated the quenching weight (3.4) for the medium-induced gluon 
energy distribution in the N = 1 opacity approximation (2.19). In general, the 
quenching weight depends in this approximation on the characteristic gluon 
energy u c , the kinematic constraint R = x 2 & c L, and the opacity n L. 
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Fig. 9. The continuous part of the quenching weight calculated in the single 
hard scattering approximation with opacity no L = 1 for a hard quark (upper row) 
or hard gluon (lower row). 



For the numerical results presented in Figs. 8 and 9, we use uqL = 1. The gluon 
energy distribution uj^ depends linearly on n L, but the quenching weight 
shows a complicated dependence on n L; it has to be calculated separately 
for each value of n L from eqs. (3.2) and (3.3). However, since n L multiplies 
the Casimir factor Cr in the gluon energy distribution, the quenching weight 
for gluons with uqL = 1 is identical to the quenching weight for quarks with 
n L = Ca/Cf = 2.25. Vice versa, the quenching weight for quarks given in 
Figs. 8 and 9 can be viewed as a quenching weight for gluons with n L = 
Cf/ Ca- 
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C. Comparison: multiple soft vs. single hard scattering approximation 

In the opacity approximation, one specifies both the average transverse mo- 
mentum squared ~ /z 2 uqL transferred to the projectile and the average num- 
ber n L of scattering centers involved in this momentum transfer. This is in 
contrast to the multiple soft scattering approximation which specifies the av- 
erage transverse momentum squared transferred to the projectile irrespective 
of the number of scattering centers involved. Thus, the single hard scattering 
approximation contains one additional model parameter, the opacity no L. 

Despite this difference, we want to compare the quenching weights obtained 
in both approximations. To this end, we start from the relations 

R ~ (n L) R , u c ~ (n L)u c , (3.12) 

discussed in section 2 C. Keeping the values of R, u c and R, Cu c fixed, we "fit" 
the opacity n L such that the quenching weights obtained in both approxima- 
tions show the best agreement. This allows to discuss for both approximations 
differences in functional shape which cannot be removed by a change of model 
parameters. 




BDMPS 
GLV (u) c =u) c ) 
GLV (w c =o) c /3) 



0.5 1 1.5 2 2.5 

AE/o; c 

Fig. 10. Comparison of the quenching weights for infinite system size 
(R, R — > oo) calculated for a hard gluon in the multiple soft ("BDMPS") and 
single hard ("GLV") scattering approximation. For rescaled characteristic gluon en- 
ergy, uj c ~ (uqL) Q c , n$L = 3, the agreement between both probability distributions 
improves, see text for further discussion. 
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We start by comparing the R- and _R-dependences of the discrete weight po 
calculated in the multiple soft (Fig. 6) and single hard (Fig. 8) scattering ap- 
proximation, respectively. For the choice R ~ 3 R, the curves show better 
agreement. However, the excess above unity for R < 100 is much more pro- 
nounced in the multiple soft scattering approximation, than the excess above 
unity for 3.R < 100 in the single hard scattering approximation. This indicates 
that the specific destructive interference effects discussed in section 3 A play 
a more important role in the multiple soft scattering approximation. 

In Fig. 10, we compare both approximations in the limit in which the constraint 
on the transverse momentum phase space is removed (i.e. R,R^ oo). For the 
opacity n L = 3, the maximum of the quenching weight takes the same value 
in both approximations. However, significant differences can be seen in the 
functional shape. The gluon energy distribution is harder in the single hard 
scattering approximation (see Fig. 5) and this is reflected in a more pronounced 
large energy tail of the quenching weight. 

We regard the remaining differences between both approximations as an indi- 
cation of the intrinsic theoretical uncertainties in evaluating the gluon energy 
distribution (2.1). 



4. Medium-induced gluon radiation for an expanding medium 

Hard partons produced in the initial stage of ultra relativistic nucleus-nucleus 
collisions are propagating through a strongly expanding medium. This re- 
sults in a time-dependence of the transport coefficient q(£) which can be 
parametrized in terms of a power law, 



The expansion parameter a determines the dynamical evolution of the 
medium: a = characterizes a static medium. A one-dimensional, boost- 
invariant longitudinal expansion is described by a = 1. This value is supported 
by hydrodynamical simulations of the early stage. In general, however, an ad- 
ditional transverse expansion can lead to larger values a < 3. The maximal 
value go of the transport coefficient is reached at the time of highest density of 
the system which is the formation time £o- This formation time may be set by 
the inverse of the saturation scale p sa t 17 , resulting in pa 0.2 fm/c at RHIC and 
pa 0.1 fm/c at LHC. The difference between £o and the production time of the 
hard parton is negligible for the calculation of the gluon energy distribution 
(2.1). It will be ignored in what follows. 




(4.1) 
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In this section, we discuss the range of validity and the form of a dynamical 
scaling law 11 which relates the gluon energy distribution (2.1) in a collision of 
arbitrary dynamical expansion to an equivalent static scenario. 



1. Multiple soft scattering in an expanding medium 

In Appendix C, we give details of the calculation of the gluon energy dis- 
tribution (2.1) for values of the expansion parameter a < 3 in the multiple 
soft scattering approximation (2.3). As reported previously 11 , we observe a 
scaling law which relates the time-dependent transport coefficient (4.1) to an 
equivalent static transport coefficient q, 




10 io to ^o^o 1 1 10 4 10 3 10 2 10 1 

0)/(ff L 2 /2) CO /*"" C0 C = CO / (ef l_ 2 / 2 ) 



Fig. 11. The gluon energy distribution calculated in the multiple soft scattering 
approximation with expansion parameter a = (solid line), a = 0.5 (dashed line), 
a = 1.0 (dotted line) and a = 1.5 (dash-dotted line). Curves in the right hand 
column are scaled according to (4.2). 

~q=^J LH °dU^-^)m)- (4-2) 
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As seen in Fig. 11, the gluon energy distributions for different values of the 
expansion parameter a differ by orders of magnitude if plotted in units of 
the same characteristic gluon energy uj c = |go L 2 and kinematic constraint 
R — uj c L. However, if plotted in units of the rescaled gluon energy \qL 2 and 
the rescaled kinematic constraint ^qL 3 , they agree approximately over a large 
parameter range. 

For practical purposes, the accuracy of the scaling law (4.2) is satisfactory 
for R > 100. Concerning the deviations from the scaling law for R < 100 
(see Fig. 11), we make the following comments: In practice, these deviations 
are negligible since p ~ 1 for R < 100 and thus no significant medium 
modification occurs. Technically, the static case (a = 0) is calculated for a 
box profile in the longitudinal density of scattering centers. On the other 
hand, in the expanding scenarios, the density profile degrades more smoothly 
with increasing path length, and the discontinuity at path length L is less 
important. The strength of destructive interferences between medium-induced 
and vacuum gluon radiation depends on this discontinuity. This may explain 
why for R = 40 the rescaled gluon energy distribution in Fig. 11 is more 
suppressed in the static case than in the expanding case. 

2. The opacity expansion for an expanding medium 

In Appendix B, we give analytical expression for the single hard scattering 
limit of the gluon energy distribution (2.1) in a medium with expansion pa- 
rameter a. The analytical form of (2.1) changes with the expansion parameter 
a. We derive an explicit expression for the Bjorken scaling case a — 1, 

7^ = 2^^ (n £o) / — Re [-Ei[-ir] + H~^] + 1e] 
duj it Jo r 

x (jl 1 

\ r + 7 ^(k 2 + r + 7 ) 2 - An 2 r 

Here, 7^ = 0.577... denotes Euler's constant and the exponential integral 
function is Ei[z] = -J^dte^/t. 

To relate the gluon energy distributions for a static medium (2.19) and a 
Bjorken scaling expansion (4.3), we determine the dynamically averaged den- 
sity of scattering centers following eq. (4.2) 

2 f L +(o 2n £ A , 

n = j2 j ^ " n (0 = — — • ( 4 - 4 ) 

This equation suggests that the gluon energy distributions in the static and 
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Bjorken expansion case show agreement if the prefactor (no£o) m eq. (4.3) is 
replaced by \nL where n determines the density of scattering centers of the 
equivalent static scenario. In Fig. 12, we test this suggestion numerically for 
different values of the kinematic constraint R. 

Remarkably, for sufficiently large kinematic constraint R > 100, the gluon en- 
ergy distribution shows dynamical scaling according to (4.4) for gluon energies 
iv > lv c . This is the dominant kinematic region in the opacity expansion (see 
the discussion in section 2B, eq. (2.21) ff). Thus, despite the deviations from 
the scaling law for tv < iv c , the logarithmically enhanced contribution to the 
average energy loss 
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Fig. 12. The gluon energy distribution calculated in the single hard scattering 
approximation for a static medium (dashed line) and for a medium with Bjorken ex- 
pansion (solid line) . Curves for the dynamically expanding case are scaled according 
to eq. (4.4). 



r dl N=1 

(AE) N - 1 = / diviv — shows scaling 18 with (4.4). 

J du 



(4.5) 



This is also known to hold in the soft multiple scattering approximation and 



22 



it is consistent with results obtained on the basis of twist-4 matrix elements 20 . 

In the region uo < Cu c , significant deviations from the scaling law (4.4) are seen 
in Fig. 12. However, the logarithmic plot overemphasizes the importance of 
these deviations. First, they occur in the sub dominant region which is less 
important for calculating the quenching weights. Second, these deviations do 
not exceed 30 percent in the physically relevant parameter range 100 < R < 
40000 in which significant medium modifications can be expected. 

For practical purposes, the scaling law (4.4) is thus satisfactory. Quenching 
weights for a dynamically expanding scenario can be obtained by calculating 
the quenching weights of the dynamically equivalent static scenario according 
to (4.4). 



5. Angular dependence of radiation probability 



The maximal angle under which a gluon can be radiated is given by the upper 
bound on the transverse momentum integral in (2.1), 

Q^^ = X . (5.1) 



Thus, for fixed values of the characteristic gluon energy u c and of the kine- 
matic constraint R = uj c L, a decreasing value of R x = x 2uJ cL determines 
the medium-induced energy radiated into a cone of opening angle 6. In 
this section, we denote explicitly the dependence of the quenching weight 
P(AE,u c , R x ) on uj c and x 2 R. This quenching weight determines the prob- 
ability that an additional energy fraction AE is radiated inside the opening 
angle x — ®- From the Figs. 1 and 3, we know that the more collinear com- 
ponent of the medium-induced gluon radiation is harder. 

For fixed values of uj c and R = u c L, the gluon energy distribution radiated 
outside the opening angle 6 is given by 

dl >e dl dl 



The probability that an additional energy fraction AE is radiated outside the 
opening angle x ~ is obtained by inserting (5.2) into the Mellin trans- 
form (3.2), (3.3). For the current work, we did not calculate this probability; 
there is no simple way to obtain it without Mellin transform directly from the 
quenching weights tabulated in section 3. 
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The calculation of the average energy loss outside an angle 6 is simpler. It 
can be calculated from the quenching weights tabulated in section 3 



r dl >@ 
(AE)(Q) = duo uo— — (uo c ,R = uo c L) 
J duo 



= dEE 



P{E,uo c ,R = uo c L)-P(E,uo c ,R x = X 2 uj c E) . (5.3) 



In Fig. 13, we compare the angular dependence of the average parton energy 
loss (5.3) in the multiple soft and single hard scattering approximation. 

In the single hard scattering approximation, the integral (5.3) diverges log- 
arithmically in the ultraviolet for % = O ^ 0. For the calculation of 
(AE) N=1 (Q = 0) in (2.22), we have cut off this divergence by limiting the 
energy radiated away to uo < E. For the plot in Fig. 13, we restrict instead the 
calculation to sufficiently large angles 6 for which the second term in (5.3) 
provides an ultraviolet cut-off. Thus, for small angles 6 < 10° where AE ~ E, 
one overestimates (AE) N=1 (Q). 

R =R=2000 




angle (degrees) angle (degrees) 



Fig. 13. The average energy loss (5.3) radiated outside an angle as calculated 
in the multiple soft (lower three lines) and single hard (upper three lines) scattering 
approximation for a jet of energy E = 100 GeV. 



In the single hard scattering approximation, the region uo > oo c is dominant. 
This hard part of the spectrum is emitted under angles 6 < ji/uo c and thus 
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appears as a logarithmic enhancement in the collinear region. In the multiple 
soft scattering approximation, however, the dominant radiative contribution 
lies in the soft region u < uj c which is emitted under relatively large opening 
angles > k±/uj c ~ q 1 ^ /oo^ A . For smaller opening angles, the average energy 
loss (AE)(Q) does not increase further. Indeed, multiple soft scattering results 
in a shift in transverse phase space which is known to deplete (AE)(Q) at very 
small angles 6 ' 21 . Thus, in the multiple soft scattering approximation, there is 
no ultraviolet divergence at small angle 6. 

To compare the single hard and multiple soft scattering approximations for 
(AE)(Q), we proceed in analogy to the discussion in section 3 C: varying 
the opacity, we find the best agreement between both approximations for 
n L = 3. A qualitative difference which cannot be adjusted by the choice of 
the additional model parameter no L persists for small angles only. Its origin 
is explained above. Thus, Fig. 13 indicates that for comparable sets of model 
parameters u c , R and Q c , R, n L respectively, the multiple soft and single 
hard scattering approximations lead to a comparable angular dependence of 
(AE)(G) for > 10°. 



O 2.5 




2 4 6 8 10 12 14 16 18 20 

(q(L/2) 3 ) 1/2 cone 

Fig. 14. Angular dependence of the average energy loss (5.3) in the multiple 
soft scattering approximation as a function of the rescaled parameter g(L/2) 3 0. 

The calculation of the angular dependence in Fig. 13 was done for quenching 
weights calculated for a static medium. In general, the dynamical scaling laws 
(4.2) and (4.4) relate these to the quenching weights of dynamically expanding 
scenarios. However, this is not the case for the small values R, R < 100 which 
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enter the calculation of the small angular dependence (0 < 10°) in eq. (5.3). 
For these small values of R, R, dynamical scaling breaks down (see Figs. 11 
and 12). Since typical jet cone openings correspond to larger angles > 10°, 
we did not make an effort to quantify the remaining dependence of (AE)(Q) 
on the collective expansion of the collision region. 

In the multiple soft scattering approximation, (AE)(Q = 0) is finite. Baier et 
al. observed 22 that the ration (AE)(Q)/ (AE)(Q = 0) is an universal quantity 
which depends solely on qL 3 Q. Fig. 14 shows that this statement remains 
approximately true in the presence of a finite kinematic constraint R. 



6. Applications of quenching weights 



In this section, we use quenching weights to calculate in two alternative ways 
the medium-induced suppression of hadronic high transverse momentum spec- 
tra. In section 6 A, we determine the quenching factor Q(p±) and in section 6 B 
we calculate medium- modified parton fragmentation functions. In section 6 C 
we finally discuss the relation of both calculations to the nuclear modification 
factor measured at RHIC. 



A. Quenching factors for hadronic spectra 



The medium-dependence of inclusive transverse momentum spectra can be 
characterized in terms of the quenching factor 1 

= to^pmi = f dAEP(AEj (n+^ifA) . (6 . 1} 

v ' da vac (p ± )/dp 2 ± J \ da vac (p±)/dp 2 ± J 



Here, the spectrum da vac (p±)/dp 2 ± is unaffected by medium effects; it is deter- 
mined e.g. in proton-proton collisions. Equation (6.1) relates it to the medium- 
modified transverse momentum spectrum da mcd (p±)/dp 2 ± measured e.g. in 
nucleus-nucleus collisions. We work in the longitudinally comoving frame in 
which the total energy of the produced parton is directed orthogonal to the 
beam. Due to the presence of the medium, a parton produced initially with 
transverse momentum p± + AE looses an additional energy AE with proba- 
bility P(AE). This defines the quenching factor (6.1). 

If one assumes a power law fall-off of the p^-spectrum, then the quenching 
factor (6.1) can be calculated explicitly, 
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Q(p ± ) ~ J dAEP(AE) 



P± 



(6.2) 



In general, the effective power n depends on p± and y/s. It is n 
kinematic range relevant for RHIC. 

q"=1 GeV 2 /fm 
n=4, L=2frn n=4, L=5fm 
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Fig. 15. The quenching factor (6.2) calculated in the soft multiple scattering 
approximation for a s = 1/2. Upper row: calculation in the R — > oo-limit but with 
a varying sharp cut-off on the infrared part of the gluon energy distribution uij^ 
determining the quenching weight. Lower row: the same calculation is insensitive to 
infrared contributions if the finite kinematic constraint R = lo c L < oo is included. 



To compare directly to published results 1 , we calculate the quenching factor 
(6.2) in Fig. 15 for parameter values used previously. The transport coefficient 
is taken to match expectations for a hot medium, q = ^-^p- ■ Given the in- 
medium path length L, this defines u c and R in the multiple soft scattering 
approximation. In the single hard scattering approximation, parameters are 
chosen for opacity n L = 1 by identifying R = R and Cu c = uj c . The effective 
power n in (6.2) is set to its asymptotic value n = 4. To be quantitatively 
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comparable with Ref. 1 , we use a s = 1/2 in Figs. 15 and 16 while all other 
numerical results are given for a s — 1/3. 

Fig. 15 shows the numerical results obtained in the multiple soft scattering 
approximation. The perturbative calculation of the gluon energy distribution 
uj^ cannot be trusted for soft gluon energies cu ~ O(Aqcd) where perturba- 
tion theory breaks down. To quantify the sensitivity of their calculation to 
this infrared region, Baier et al. 1 introduced a sharp cut-off on the R — > oo 
gluon energy distribution which was varied between u cut = and u cut = 500 
MeV. The resulting uncertainty is seen as a ~ 20% variation of the quench- 
ing factor Q(p±) in the upper row of Fig. 15. The finite kinematic constraint 
R = uj c L < oo depletes the infrared region of the medium-induced gluon ra- 
diation spectrum (see Fig. 1). This constraint is a generic consequence of a 
finite in-medium path length; it removes almost completely the sensitivity of 
the calculation to the uncontrolled infrared region (see Fig. 15, lower row). Re- 
markably, it also tends to flatten the p^-dependence of the quenching factor. 
We shall return to the consequences of this observation when we zoom into 
the region p± < 10 GeV in section 6 C. 
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Fig. 16. The same quenching factor (6.2) as in Fig. 15, here calculated in the 
single hard scattering approximation. 
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We have calculated the quenching factor (6.1) in the single hard scattering 
approximation, see Fig. 16. Since the dominant contribution comes in this case 
from the hard part of the spectrum, u > uj c , the sensitivity to the infrared 
cut uj cut is much reduced in comparison to Fig. 15. However, realistic finite 
kinematic constraints R remove a much larger part of the soft spectrum. As 
a consequence, the absolute value of Q(p±) increases significantly if finite 
kinematic constraints are imposed, and the p_i_-dependence tends to flatten. 

In section 3 C, we observed that for quenching weights the best agreement 
between single hard and multiple soft scattering approximation is for u c = 3uj c , 
R = 3R, n L = 3. Remaining differences come from the fact that the single 
hard scattering approximation shows a dominant contribution for u > oo c while 
the multiple soft scattering approximation shows a dominant contribution for 
uo < uo c . However, the quenching factor Q(p±) for small p± is sensitive only to 
the soft region uj < uj c , (uj < uj c ) in both approximations. This is so because 
medium-induced gluons cannot carry away more than the total energy E q 
of the parent parton, and hence uj < E q < uj c (cu c ) at small p±. Thus, the 
simple relation uj c = 3uj c , R = 3R does not hold for p± < uo c . This is seen 
in Fig. 17. For small p±, the multiple soft scattering approximation results in 
a much stronger suppression than the single hard one calculated for rescaled 
parameters u c = 3u c , R = 3R. 
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Fig. 17. Comparison of the quenching factor (6.2) calculated for a s = 1/3 and 
effective power n = 7 in the multiple soft, single hard and rescaled single hard 
scattering approximation. 
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B. Medium-modified fragmentation functions 



Medium-induced gluon radiation off hard partons modifies the fragmentation 
and hadronization of final state partons, thus affecting hadronic p_i_-spectra. In 
section 6 A, we calculated this effect in terms of the quenching factor Q(p±). 
Alternatively, this quenching factor can be determined from medium-modified 
fragmentation functions, which we discuss now. 

In the Q CD-improved parton model, hadronic cross sections for high-p^ 
hadropro duct ion are calculated by convoluting the perturbatively calculable 
hard partonic cross section da q and the (final state) fragmentation function 



For notational simplicity, we do not denote the additional convolution of do~ q 
with the (initial state) parton distributions. The leading hadron h carries an 
energy fraction z = Eh/ Q of the total virtuality of the partonic collision, which 
is a fraction x = Eh/E q of the energy of its parent parton. The parent parton 
carries the energy fraction y = E q /Q. 

If the parent parton loses with probability P(e) an additional energy frac- 
tion e = %^ prior to hadronization, then the leading hadron is a fragment 
of a parton with lower energy (1 — e)E q ; thus, it carries a larger fraction 
of the initial parton energy. The inclusion of this effect amounts to replac- 
ing the fragmentation function D h / q (x,Q 2 ) in (6.3) by the medium-modified 
fragmentation function 23 ' 24 



To calculate eq. (6.4), we use the recent LO fragmentation functions of 
Kniehl, Kramer and Potter 25 ("KKP"). These improve over previously avail- 
able parametrizations 26 . However, the KKP parametrization still shows signif- 
icant uncertainties in the large- a; region relevant for hadronic p^-spectra 27 . For 
alternative approaches towards medium-modified fragmentation functions, see 



D h / q (x,Q 2 ), 




(6.3) 



(6.4) 



Refs. 28 ' 29 . 
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Fig. 18. The medium-modified fragmentation function (6.4) for Q 2 = (10 GeV) 2 
calculated in the multiple soft scattering approximation for a medium of length 
L = 6 fm. 



We have calculated the medium-modified q — > 7r and g — > 7r fragmentation 
functions (6.4), using the quenching weights in the multiple soft (Fig. 18) and 
single hard (Fig. 19) scattering approximation. The energy of the parent parton 
is set to the virtuality of the hard process 11 , E q ~ Q. The medium-induced 
fragmentation functions decrease with increasing density of the medium since 
the probability of a parton of initial energy E q to fragment into a hadron of 
large energy xE q decreases with increasing parton energy loss. They should 
be trusted for sufficiently large momentum fractions (x > 0.1 say) only. The 
reason is that the hadronized remnants of the medium-induced soft radiation 
are not included in the definition of (6.4). These remnants are soft - they can be 
expected to give an additional contribution in the region x < 0.1. The neglect 
of these remnants in (6.4) implies that the normalization of D^ d \x, Q 2 ) is a 
factor / deeP(e) too small, 

J* dxxD^\x) ~ J* dxxD h/q (x) J de (1 - e) P(e) . (6.5) 

For the suppression of high-p^ hadronic spectra, this normalization error is 
unimportant since the main contribution comes from the region of larger x. 
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To estimate the corresponding suppression of high-p^ hadronic spectra, we ex- 
ploit that the fragmentation function in (6.3) is weighed by the partonic cross 
section da q /dpj_ (we work in the frame in which the total energy of the pro- 
duced parton goes transverse to the beam). In the kinematic range relevant for 
RHIC (v/i = 200 GeV and P± ~ 10 GeV), one finds 30 d<f/dp[ ~ l/p^" PJ - ) 
with n(y/s,p ± ) ~ 6. Thus, eq. (6.3) effectively tests x n ^' p ^D { ™ d) (x,Q 2 ). 
The suppression factor 



R ff(P±) = 



6 n (med)/ 
•^max^h/q l^max, 



X max^h/q( x max: P±) 



(6.6) 



P_L=X n 



provides a simple estimate of the reduction of hadronic p_i_-spectra. In eq. 
(6.6), x max denotes the maximum of x n ^' p± ^D^ d \x, Q 2 ) and corresponds 
to the most likely energy fraction p± = x max E g of the leading hadron. The 
suppression factor can be read off easily from the lower rows of Figs. 18 and 
19. We now compare this suppression factor to the quenching factor Q(p±) in 
(6.2). 
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C. The nuclear modification factor 



Experimental situation: Published data for Au+Au collisions at v /snn = 130 
GeV show for p± < 6 GeV a suppression of neutral pion 31 and charged 
hadron 31 ' 32 transverse momentum spectra if compared to spectra in p+p col- 
lisions rescaled by the number of binary collisions. This suppression is most 
pronounced (up to a factor ~ 5) in central Au+Au collisions and smoothly 
approaches the binary scaling case with decreasing centrality. Within error 
bars, the suppression factors of 7r° and charged hadron spectra agree, though 
central values for the suppression of 7r° production are slightly lower 31 . In ad- 
dition, a maximal azimuthal anisotropy ^(rl) of hadroproduction is found to 
persist up to the highest transverse moment um 33 ~ 35 . These data indicate the 
importance of final state medium effects up to p± < 6 GeV. 

Preliminary data shown at the Quark Matter 2002 conference confirm these 
findings for Au+Au collisions at ^snn = 200 GeV; they extend many observa- 
tions up to p± ~ 10 GeV. In particular, data for the nuclear modification fac- 
tor show an approximately constant maximal suppression within 6 < p± < 12 
GeV for charged hadrons 36-38 and up to p± < 8 GeV for n° spectra 39 ' 40 . The 
azimuthal anisotropy V2(p±) of charged hadrons stays close to maximal up to 
p± < 10 GeV 41 . Moreover, the disappearance of back-to-back high-pj_ hadron 
correlations 42-44 provides an additional indication that final state medium ef- 
fects play a decisive role in hadroproduction up to p± ~ 10 GeV. 

Theoretical situation: Parton energy loss has been proposed to account for the 
small nuclear modification factor 16 ' 20 ' 45 , the azimuthal anisotropy 47 ' 18 ' 48 ' 46 and 
the disappearance of dijets 50 ' 51 . Quantitative studies indicate, however, that 
in the kinematic regime relevant for RHIC (p± < 12 GeV), p_|_-broadening 51 ' 20 , 
shadowing 45 ' 20 , formation time 16 and possibly other effects contribute to the 
high-px nuclear modification as well. Models have been proposed which ac- 
count for hadronic quenching without taking recourse to parton energy loss. 
Instead, these models invoke string percolation 52 , small hadronization time 
arguments 53 , saturation physics 54 , the dominance of parton recombination 
over parton fragmentation 55 , or initial state formation time arguments 56 . The 
consistency and applicability of these models is currently under debate. For 
the p_i_-range accessible to RHIC, the competing hadronic effects may make 
it difficult to disentangle quantitatively the contribution of parton energy loss 
from the measured hadronic suppression pattern. The transverse phase space 
accessible to LHC (E± < 200 GeV) may turn out to be a qualitative advantage 
with this respect. 

Model comparison: In Fig. 20, we compare the two definitions (6.2) and (6.6) 
of the quenching factor to the nuclear modification factor measured by the 
PHENIX Collaboration in the 7r°-spectra 39 ' 40 of central Au+Au collisions at 
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^/snn = 200 GeV. We do not include the nuclear modification factor for 
charged hadrons 36-38 in Fig. 20, since charged hadrons are likely to be domi- 
nated at high p± by baryons whose production mechanism may involve addi- 
tional non-perturbative effects 57 . 
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Fig. 20. The nuclear modification factor for 7r°-production 39 ' 40 compared to 
model calculations involving parton energy loss only. The lines show the quenching 
factor (6.2) with effective power n = 7, and the suppression factor (6.6) calculated 
from medium-modified fragmentation functions. They are given in the limiting cases 
where all parent partons are either quarks (upper lines) or gluons (lower lines). 
Calculations in the multiple soft scattering approximation use R = 2000, uj c = 67.5 
GeV, corresponding to q = 0.75 G ^ and L = 6 fm. In the single hard scattering 
approximation, we use R = R, cd c = lo c . 



Our comparison makes several simplifying assumptions which can be improved 
in further studies: i) We do not model a realistic space-time geometry of the 
collision. Instead, we work for a fixed in-medium path length L = 6 fm. As a 
consequence, we do not compare to the centrality dependence of the nuclear 
modification factor for which a realistic distribution of in-medium path lengths 
and their impact parameter dependence is needed, ii) We do not calculate the 
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hard partonic matrix elements entering e.g. (6.3). Instead, we determine the 
quenching factors directly from (6.2) and (6.6). As a consequence, we do not 
know the p_i_-dependent fractions of hard quarks and of hard gluons which 
fragment into a leading ir°. Since parton energy loss is different for quarks and 
gluons, we plot in Fig. 20 the limits for which all parent partons are quarks or 
gluons, respectively. The realistic curve lies in between these limits. Since the 
ratio of parent quarks over parent gluons increases with p±, the p_i_-dependence 
of the realistic curve will be slightly steeper than the limiting cases presented 
in Fig. 20. 

From Fig.20, we can draw several conclusions. First, the two definitions of 
quenching factors in terms of hadronic spectra (6.2) and in terms of fragmen- 
tation functions (6.6) lead to quantitatively comparable results. Second, the 
multiple soft and single hard scattering approximations for parton energy loss 
lead to quantitatively comparable results. The slight variations in pj_-slope 
should be regarded as theoretical uncertainties in approximating (2.1) and 
do not give any preference to either approximation scheme. Third, a calcula- 
tion based on partonic energy loss only can reproduce the magnitude of the 
observed nuclear modification factor. Moreover, it results in a very shallow in- 
dependence which seems consistent with the current quality of experimental 
data. For an interpretation of the model parameters used in Fig.20, we use for 
a Bjorken scaling expansion (a = 1) the relation 11 ' 18 between the kinematic 
constraint R and the initially produced gluon density, 



L 2 dN9 
R\^ 



R=T»—> ( 6 - 7 ) 



where Ra denotes the nuclear radius. The extracted value ~ 2000 is ap- 
proximately a factor 2 larger than a previous estimate 18 based on data from 
elliptic flow. Given the theoretical uncertainties of parton energy loss calcula- 
tions below p± < 10 GeV, this factor 2 mismatch constitutes no inconsistency. 



7. Conclusion 



In recent years, phenomenological applications of medium-induced parton en- 
ergy loss were based mainly on two different approximations of the medium- 
induced gluon energy distribution: the multiple soft BDMPS-Z scattering ap- 
proximation and the opacity approximation. It remained unclear, however, 
to what extent these approximations differ. Also, the quenching weights were 
studied mainly in kinematic limiting cases. Here, we have studied in both 
approximations the medium-induced energy distributions (section 2), the cor- 
responding quenching weights (section 3), and the extension of these calcula- 
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tions to dynamically expanding collision region (section 4) and to the angular 
dependence of the medium- induced radiation pattern (section 5). The single 
hard scattering approximation is dominated by the hard region uo > uo c of the 
gluon energy distribution while the multiple soft scattering approximation is 
dominated by the soft region uj < u c . Despite this difference, both approxi- 
mations lead to quantitatively comparable results if comparable sets of model 
parameters are used. Numerically, we determine the correspondence 

w c ~3w c , R~3R, (7.1) 

which relates the BDMPS transport coefficient q and the Debye screening 
mass u 2 via an opacity no L ~ 3, see eq. (3.12). Deviations from (7.1) can be 
understood in terms of kinematic constraints on the hard part of the gluon 
energy distribution (see discussion of Fig. 17). 

The main result of this paper is the calculation of quenching weights in section 
3. We explained how to calculate from these quenching weights the nuclear 
modification of hadronic spectra. There are indications that the interpretation 
of RHIC data on hadronic quenching requires additional physics effects beyond 
the parton energy loss (see discussion in section 6 C). However, to discriminate 
energy loss contributions from these additional effects, a quantitatively reli- 
able discussion of the current theoretical calculations of parton energy loss is 
needed. We hope that the CPU-inexpensive subroutine for quenching weights 
which accompanies this paper will prove a valuable tool to this end. Also, this 
routine can be used to explore observable consequences of parton energy loss 
in nucleus-nucleus collisions at LHC or for processes in cold nuclear matter. 
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A. The BDMPS-limit of the energy distribution (2.1) 

Here, we establish that the R — > oo limit of the gluon energy distribution 
(2.1) coincides with the BDMPS result (2.9). Using the saddle point approx- 
imation (2.3), the energy distribution (2.1) can be written in the form given 
in equations (All) and (A12) of Ref. 6 . Integrated over transverse momentum 
< k± < x^, one finds 
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CO 



dl ct s Cp 

duJ 7T 
h 



(h + h) , 



^2Re£ d yi J^ dm 16 A 
iA 4 B 4 { X uo) 2 \ 



4(D 4 - iA 4 B 4 ) 2 



1 r L 
-Re / dy t 

CO Jo 



4A_ 5 

B R 



cxp 



exp -j 
AA 5 B 5 



A(D 4 - iA 4 B 4 ) 2 
-1 , 



- 1 



(A.l) 

, (A.2) 
(A.3) 



where 
A 4 
D 4 
A, 



coQ 



2sin(Q(y L -y L )) ' 
-n C(L - y L ) , 



B 4 = cos(Q(y L -y L )) , 



uoQ 



2sin(fh/ L ) 



B 5 = cos(fh/ L ) , 



and 



w 4a; 



(A.4) 
(A.5) 
(A.6) 



(A.7) 



The limit R — > oo is obtained by taking x — > oo in eqs. (A.2) and (A.3): 
lim I 4 = —2Re / dyi / cft/j 

-R— >oo jo Jt/, 



sin (ft (y, - j/,)) 
- ft 2 



JO JO SI 

= 2Re / d yi fl- 
Jo si 



sin 2 (fi (yi -yi)) 



R 



COS (ffy;) 

sin (Qyi) ' 

lim I 5 = —2Re / dy^- 



sin (fh/ ; ) cos (ft^ 



(A.8) 
(A.9) 



Both integrals are logarithmically divergent but this divergence cancels in the 
sum 



lim (I 4 + h) = 2Re In [cos (ftL)] . 



(A.10) 



This coincides with the BDMPS result (2.9). 
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B. Gluon energy distribution to first order in opacity 



In this appendix, we calculate the first order in opacity of the gluon energy dis- 
tribution (2.1) for a Yukawa-type elastic scattering center with Debye screen- 
ing mass jj,: 



Kq)l 2 = 



7r(q 2 + n 2 ) 2 



(B.l) 



According to eq. (6.4) of Ref. 4 , the energy distribution takes the form: 



dl 



N=l 



00- 



a, 2C R fx" 



du (2vr) 2 u 2 



r dk r dq / k • q Z(Q, , (B.2) 

Jo Jo n(q 2 + fi 2 ) 2 



where 



k 2 

« = !;• 



Qi 



2uo 



(B.3) 



and 



Z(Q,Q 1 ) = limRe f° dy f 



-ey-ey 



y 



(B.4) 
(B.5) 



Irrespective of the value of the expansion parameter a in (B.4), the expression 
factorizes in the form (B.5). In order to simplify (B.2), we shift the integration 

variables by q — > q — k, k — > k-J^f and q — > qWy. This leads to 



dI N = l A a s C R r« f 2 * 
uj — = 4 



PK r Z7T POO Q 

dk dip qdq — 
Jo Jo Jo ok 



du 



71 



k 2 + 2 kq cos p + q 2 + 7 / 



(B.6) 



where 



7 = 



H 2 L 

2oo ' 



R_ 

2l 



(B.7) 
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The ip- and /c-integration in (B.6) can be done analytically, 
u— = 4 -/ qdqZ(q / L) 



duj n L Jo 



x I 7 — - -r— ) • (B.8) 

' '(k 2 + q 2 + 7 ) 2 - 4/tV 9 2 + 7/ 



1. For the static case, a — 0, the phase factor (B.5) reads 
-LQi + sin(LQi) 



z a =o(Qi) — 



(B.9) 



and 



a;- 



,« S C R roc 

=4 n L 7 / gdg 

da; 7T jo 



00 q 2 - sin(q 2 ) 



x 



/ + 7 ^J(k 2 + q 2 + 7 ) 2 - 4k 2 ? 2 



(B.10) 



Substituting r = g 2 , we find (2.19). 
2. For the Bjorken scaling case, a — 1, the phase (B.5) reads 



^=i(Qi) = 7f Re 

VI 



(Ei[-iQx(L + 6)] - Eif-ig^o]) + In 



£ + 6) 



Re [Eif-iQiL] - lnf-iQxL] - 7 £ + O(W^)] 



= a Re 



X 



iQiL e ~* - 1 

<**_-_ + O(£ /L) 



(B.H) 



Here, 7^ = 0.577. . . is the Euler constant and the exponential integral 
function Ei is defined in the text following eq. (4.3). Corrections of order 
O(£o/L) can be ignored since the time of production £ is much smaller 
than the in-medium path length L. With this approximation, one has 



u — ?=±- = 4 



du 



Ur ( c \ f 
K£o) / 

7T JO 



qdq 



Rc 



x 



7 



-Ei[-iq 2 ] +ln[-iq 2 }+-f E 

7 



q 2 + 7 



(k 2 + g 2 + 7) 2 — 4k 2 o 2 



B.12) 



Substituting r = g 2 , we find (4.3). 
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C. The dipole approximation for an expanding medium 



In this appendix, we follow Ref. 19 in giving explicit expressions for the path 
integral (2.14) in the dipole approximation, 



Drexp i- J d£ 



,2 



r — 



0) r 2 



(C.l) 



Equation (C.l) is the path integral of a 2-dimensional harmonic oscillator with 
time-dependent imaginary frequency 



(C.2) 



and mass uj. The solution of (C.l) can be written in the form 
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£(ri,j/i; r 2 ,y 2 |u;) 



2iriD(y u y 2 ) 



exp [z5ci(ri,2/i;r2,2/2)] 



(C.3) 



Here, the classical action S c \ in (C.3) takes the form 



5 c i(ri,yi;r 2 ,y 2 ) = - 



r c i(0 • |r cl (0 



V2 



(C.4) 



where the classical path r c i(£) satisfies the homogeneous differential equation 



d 2 mo) 



r c i(0 = 0, 



(C.5) 



with initial conditions 

r c i(yi) = ri , and r d (y 2 ) = r 2 . 



(C.6) 



The fluctuation determinant D in (C.3) satisfies 

d 2 mo) 



de 



(C.7) 



with initial conditions 



d 



D(e,0 = 0, and -D(t,0\z=?= 1- 



(Ci 
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In practice, D is found by combining the two independent (scalar) solu- 
tions /i, f 2 of (C.5), 



and fixing the norm M by the initial condition (C.8). The solution of (C.l) 
can be written in terms of £>(£,£') and two £- and ^'-dependent variables ci, 

C2, 



We consider three cases: 

1. The case a < 2: 

For this case, explicit expressions for (CIO) are given in Appendix B 
of Ref. 19 . The two independent solutions of the homogeneous differential 
equation (C.5) are 



tf(6O=^(/i(O/2(O-/2(0/i(O), 



(C.9) 



JC(r 1 ,y 1 ;r 2 ,y 2 \uj) 



2irD(y 1 ,y 2 ) 





(C.11) 
(C.12) 



where I v and K v are modified Bessel functions with argument 



1 



(C.13) 



2 - a 



In terms of the variable [use f2o.(£o) 




z(0 = 2uQ a (^ , 



(C.14) 



the solution (CIO) is given by 19 [we use z = z(£), z' = z(£')] 



= m nlc ^ (**7 \h{z)K v {z') - K„(z)I u (z')} , 




(C.15) 



(C.16) 



(C.17) 
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2. The cae 



ise a 



2: 



In this case, the two independent solutions of the homogeneous differ- 
ential equation (C.5) are /i(£) = £2( 1-j4 ) and /2(C) — ^( 1+A ) where 



A = yjl + 4fi2 =2 (^ ). From this, one finds 



Cl 



c 2 



1 + A / £ \ 5^ A ) ]_ — A 



(c.i8) 

(C.I9) 
(C.20) 



3. The case a > 2: 



In this case, the solution (C.I5) - (C.17) has the argument 

z(t) = 2\v\n a (s )& . 



(C.21) 



Modified Bessel functions with negative index can be avoided with 
the help of the identities: K v (z) = K-„(z) and I v {z) — I- V {z) = 
-2 s -^K u (z). 



£(£,£') = 2vy/Z? [J^)^') - , 



(C.22) 



c 1 = 2\v\^&n a (Z )[I v _ 1 (z)K v W+K v - 1 (z)I v W] , (C.23) 
c 2 = 2\u\^^ ^a(Co) + . (C.24) 



With the solution (CIO), the radiation spectrum can be written as the sum 
of three contributions 6 



<*/<**> _ Ida _a +n- dli 
—, — — —3— — —zL>f\1a + h + 1%) — — 3— 
duo 00 duo 7T duo 



vu; dl 
+ — 



da; 



(C.25) 



Here, 1$ = l/k\ is the medium-independent vacuum gluon energy distribution. 

dJ_ 

dw 



I4 and I5 determine the medium-induced part jf- studied in this paper. They 



can be computed in the dipole approximation: 



+ 



iA\B A k 2 > 
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where 



A 4 = , B 4 = Cl {y hyi ) , D 4 = - d£n(Ocj(r) (C.28) 



In the case a — 0, the functions I±i/2(z) and ^±1/2(2) entering (C.10) have 
explicit expressions in terms of exponentials. One recovers the known expres- 
sions for the static case 3 ' 6 from 

B 4 = cos(Q(yi - yi)) , Z) 4 = Jn C(L - y t ) 



2sm(Q(y l -y l )) 2 
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